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The joint effect of the density dependence of the symmetry energy and strangeness content on the 
structure of cold neutron stars is studied within the framework of a relativistic mean field theory. It 
is shown that 2Mq are obtained for repulsive YY interaction and preferably for a small or a large 
slope L. An attractive £ potential in nuclear matter will favor the appearance of strangeness in 
stars with a mass as small as ~ lMg. if, however it is repulsive only stars with a mass > IAMq 
will contain strangeness. The joint effect of reducing the symmetry energy slope and including 
hyperons is to farther reduce the radius. Neutron star maximum mass evolve non-monotonically 
with the symmetry energy slope, and the smallest masses are obtained for values L ~ 80 MeV. 
Other neutron star variables evolve nonlinearly with the slope of the symmetry energy and depend 
on the hyperon-nucleon and hyperon-hyperon couplings. The radius of a neutron star is linearly 
correlated with the neutron star total strangeness fraction and the slope is independent of the slope 
of the symmetry energy and the mass of the star. 

PACS numbers: 21.65.Ef,26.60.-c,97.60. Jd 


I. INTRODUCTION 

The structure of neutron stars depends strongly on the 
equation of state (EOS) of nuclear matter at suprasatu- 
ration densities [l] . The densities in the center of neutron 
stars can go up to 8 times the nuclear saturation density. 
Consequently, the nucleon Fermi energy rises above the 
rest mass of hyperons thus making their appearance in 
the inner layers of the star possible. Hyperons were first 
taken into account in the description of stellar matter in 
[2]. Within a relativistic mean field (RMF) approach, 
hyperons have been first included in the EOS of stel¬ 
lar matter in n®. In [H besides the usual non-strange 
mesons a , w and p , the vector meson cj) and the scalar 
meson a* with hidden strangeness were also included. 

The recent observations of the high mass neutron 
stars PSR J1614-2230 0 with 1.97 ± 0.04 M 0 and PSR 
J0348+0432 [§j with 2.01 ± 0.04 M 0 raises the ques¬ 
tion whether the interior of compact star contains ex¬ 
otic degrees of freedom, in particular, hyperons, kaon 
condensates or quark matter. Calculations using the mi¬ 
croscopic non-relativistic Brueckner-Hartree-Fock (BHF) 
formalism taking into account free and interacting hy¬ 
perons show that hyperons greatly soften the EOS, sig¬ 
nificantly reducing neutron star masses, barely reaching 
the ‘canonical’ 1.4— 1.5M 0 neutron star mass|9]]. The 
authors of 0, making use of reasonable assumptions, 
and complementing the BHF formalism with a density 
dependent Skyrme-like term mimicking many-body in¬ 
teractions Ef to test the effect of three-body hyperonic 
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forces, were only able to reproduce up to 1.6 M 0 strange 
stars. This mass will go up to 1.7 tf n if a repulsive 
interaction is taken for the E hyperon jl2| . However, 
in recent calculations u sing an auxiliary field diffusion 
Monte Carlo calculation [13J , the authors have found that 
a strong repulsive three-body force, is needed to realis¬ 
tically describe the separation energy of the A-hyperon 
from hypernuclei, within a non-relativistic Hamiltonian. 
This repulsive ANN force produces a stiff enough equa¬ 
tion of state of hyperneutron matter that satisfies the 2 
Mq constraint, shifting the onset of hyperons to densi¬ 
ties above 0.56 fm -3 [14|. Other studies were developed 
within a relativistic mean-field (RMF) approach to dis¬ 
cuss whether hyperons or other deg rees of freedom are 
possible inside neutron stars llbl|23I | . Masses of ~ 2M 0 
were only achieved by including the hidden-strangeness 
vector-meson (j>. 

The symmetry energy is relatively well constrained at 
nuclear saturation density [24[. However, an accurate 
characterization of this property at all densities is needed 
to properly describe asymmetric nuclear matter and stel¬ 
lar matter. Neutron stars are very neutron rich so we 
can expect that the symmetry energy plays a significant 
role in this context. Many efforts have been made to 
understand the density dependence of the symmetry en¬ 
ergy (e sym ), see for instance [25l - l28| . Nuclear models 
are fitted to nuclear properties and this imposes cor¬ 
relations between several nuclear properties at satura¬ 
tion density, such as the symmetry energy at saturation 
J = £sym(/Oo); the slope of the symmetry energy at satura¬ 
tion L , the curvature of the symmetry energy K sym , and 
K t a term which characterizes the isospin dependence 
of the incompressibility at saturation and subsaturation 
densities, HH- In the present work we investigate 

how the interplay between the density dependence of the 
symmetry energy and the hyperon-nucleon and hyperon- 
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hyperon interactions define the structure and strangeness 
content of a neutron star. We will work in the frame¬ 
work of a relativistic mean field approach. The density 
dependence of the symmetry energy will be modelized 
by introducing in the lagrangian density nonlinear terms 
that include the vector isovector meson p, in particular, 
we will include pc r and pu> non-linear terms [30h32[ and 
will allow the slope of the symmetry energy at saturation 
to change from 50 to 110 MeV. This procedure automat¬ 
ically establishes relations within the model between the 
different isovector properties that will be discussed. 

The influence of density dependence of the symmetry 
energy on the strangeness, central density, gravitational 
mass and radius of neutron stars will then be studied. 
Due to the lack of information on the hyperon-hyperon 
and hyperon-nuclcon interaction we will consider a set of 
hyperonic parametrizations which will take into account 
the existing uncertainties. We will consider parametriza¬ 
tions which take into account YY interaction through the 
inclusion of mesons with hidden strangeness a* and (f> []. 
This will allow to include extra repulsion at large densi¬ 
ties if the coupling of the <7* meson to hyperons is not 
too strong. 

In section EU we present the formalism, in section m 
we show the results, in section ns we draw some conclu¬ 
sions, and in appendix tables with several neutron star 
properties are presented. 

II. FORMALISM 

In the present section we present a brief review of the 
model and discuss the choice of the parameters. 


A. EQUATION OF STATE 

For the description of the EOS of neutron star matter, 
we adopt an RMF approach in which the nuclear interac¬ 
tion is described by the exchange of mesons. The baryons 
considered in this work are nucleons (n and p) and hyper¬ 
ons (A, E, and S). The exchanged mesons include scalar 
and vector mesons (a and w), isovector meson (p), and 
two additional strangeness mesons, the scalar a* and the 
vector </>. The Lagrangian density includes several nonlin¬ 
ear terms in order to describe adequately the saturation 
and high density properties of nuclear matter. In the 
present work, we consider a Lagrangian with nonlinear 
oj-p and <7-p couplings, characterized by the two coupling 
constants and A CT respectively, as introduced in [301 - 
[32| . These couplings allow us to study the effect of the 
density dependence of the symmetry energy. 

Heavy ion collision experiments [33] suggest that the 
high density EOS should not be too stiff. In the RMF 
approach these constraints on the high density EOS may 
be implemented including a fourth order term on the vec¬ 
tor isoscalar meson ui in the Lagrangian de nsity pul . l35i] , 
or including density dependent couplings [361 ]. In the 
present study we include a fourth order w term and will 
impose that at high baryonic density the EOS satisfies 
the constraints proposed in [33j, even though theses con¬ 
straints should be taken with care since they carry some 
uncertainties due to the modelling of heavy ion flow. 

For neutron star matter consisting of neutral mixture 
of baryons and leptons in /3 equilibrium, we start from 
the effective Lagrangian density of the nonlinear Walecka 
model (NLWM) 
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and m* B = m B — g<r B & ~ 9°s' 7 * the baryon effec¬ 

tive mass, gi B are the i = a,uj,p meson coupling con¬ 
stants, and *1 'i are the Dirac fields for the baryons 
and the leptons, respectively. The baryon mass and 
the lepton mass are denoted by ms and mi, respec¬ 
tively. The sum in B is over the eight lightest baryons 
n,p, A, E + , E°, E _ , S°, 5 - . The sum in l is over the two 
leptons ( e~ and g~). The constants k and A are the 
couplings of the scalar self-interaction terms and t b is 
the isospin operator. The operators £l pl , = <9 M u;„ — d v u} p , 
fiu — d p (f)ii d v (\> p and — d p p v d v p p g p [p p x /9^) 


are the mesonic field tensors. The hidden-strangeness 
mesons are represented by the <7* and the </> M fields. In 
the RMF approximation, the meson fields are replaced by 
their expectation values in the groundstate of the system 
under consideration. The Eulcr-Lagrange equations ap¬ 
plied to Eq. m and using the mean-field approximation, 
yield the meson field equations written as 





(2a) 






3 


Wo = 

\ ' 9^b ( pB\3 
to* 2 3tt 2 ' F 

(2b) 

P03 = 


(2c) 


P B 


_* 

a 0 — 

1 f kF m * B k 2 

(2d) 

m 2 a , p io V^ 2 + to;) 2 

4>o = 


(2e) 


where kp is the B baryon Fermi momentum, and the 
meson ’’effective” masses (m* 2 , i = a,LU,p) are defined 
as 


*2 

m CT 

= ml + ~u 0 + - 2A a g 2 a g 2 p plz, 
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The extra nonlinear term a-p , in the Lagrangian density, 
causes a decrease of the effective mass of the cr-meson 
with density and the uj-p term causes an increase of the 
ui meson effective mass. However, both a-p and ui-p non¬ 
linear terms increase the p meson ’’effective mass” with 
density. Consequently, the a meson field will harden at 


large densities whilst the ui and p fields will soften at 
large densities. 

For neutron star matter composed by a neutral mix¬ 
ture of baryons and leptons, we force the generalized /? 
equilibrium with no neutrino-trapping and the charge 
neutrality which can be written as following 

Pi~b z p n = qipi , p s = 0 (4a) 

Pp+ + PE+ = Pe- + P\i- + PE~ + Pa~ (4b) 

where q t is the electrical charge, bi is the baryon number, 
Pi and pi = (kp) 3 /3ir 2 are the chemical potential and the 
number density of species i. respectively. The chemical 
potentials for the baryons (pb) and for the leptons (p{) 
are defined by 

Pb = \j (kf) 2 + m\ j 2 + g UB u} o + ^y-T 3 p 0 3 

+ 50b <A), (5a) 

Pi = \J {kp) 2 + m 2 . (5b) 

The coupled (HQ), (l2bli . ©, G3K G3b mil), SB equa¬ 
tions are solved self-consistently at a given baryon density 
pB = p P + Pn + PA + Pe+ + Pe° + Pe- + P-° + Pe- • From 
the energy-momentum tensor, we calculate the energy 
density (e) and the pressure (P): 



is the mass inside radius r. The total strangeness content 
of a given star, Ng is given by 

f R r 2 

N s = 4tt / ps dr (9) 

Jo ! X M(r) 

where 

Ps = Pa + Pe+ + Pe° + Pe- + 2 psp + 2 p s - (10) 
is the total strangeness density. 


where k l F is the l lepton Fermi momentum. To calcu¬ 
late neutron star structure we solve the general relativity 
Tolmann-Oppenheimer-Volkov (TOV) pair of equations 
which are written as 


= - qr- 2M(r)] < M(r) + (8a) 

M(r)=4n f e(r)r 2 dr (8b) 

Jo 

where e is the energy density, P is the pressure and M (r) 
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B. MODEL PARAMETERS 

The RMF models are generally characterized by a set 
of nuclear matter properties at saturation density, includ¬ 
ing the nuclear saturation density po, the binding energy 
per baryon number B/A, the nucleon effective mass m*, 
the nuclear incompressibility Kq, the symmetry energy 
J and its slope L. In this work we use the TM1-2 pa¬ 
rameter set I20tl, a modified parametrization of the TM1 
parameter set [341 ] . and whose saturation properties are 
listed in table Q] At suprasaturation densities the TM1- 
2 model is slightly stiffer than the TM1 model but still 
overlaps heavy ion flow data [33j . 


TABLE I: Saturation properties of the TM1-2 model. All 
units are in MeV. 


Po (fm 3 ) 

B/A 

m* 

K 0 

Qo 

0.145 

-16.38 

595.63 

277.03 

-197.48 

Csym 

L 

■^sym 

Q sym 

k t 

36.84 

111.27 

41.90 

-32.72 

-546.27 


The parameters of the model are the nucleon mass 
niN = 938 MeV, the masses of mesons m a = 
511.198MeV, ro w = 783 MeV, m p = 770 MeV, and the 
coupling constants, which are listed in Tables HT1 and m 


The range of values covered by J and L are, respectively, 
31.32 < J < 36.84 MeV and 50 < L < 111 MeV, see 
Fig. [2] (a). These values cover the values obtained from 
isospin diffusion observables in heavy ion collisions [39} or 
the mean N/Z distributions of emitted fragments with 
radioactive ion beams [40], but are not as low as the ones 
obtained from chiral effective field theory [4lj . 



FIG. 1: (Color online) Symmetry energy versus density for 
the TM1-2 parameter set and comparison between the ujp and 
up behaviors for L = 75 MeV (thick colored lines) and L = 50 
MeV (thin colored lines). 


TABLE II: Coupling constants and meson masses for the 
TM1-2 model. 


V m a > 

f Qu, \2 

(^) 2 

A 

■5 

fm 2 

fm 2 

fm 2 



14.8942 

9.9285 

5.6363 3.52353 

-47.36246 

0.01167 


1. Density dependence of symmetry energy 

The cop and ap nonlinear terms affect the density de¬ 
pendence of the symmetry energy [301, [H, H?}, the sym¬ 
metry energy being given by 


6 yjk 2 p + ( m N - g a a 0 ) 2 

+ ——r, --- ( 11 ) 

^ rn 2 p + 2g 2 p [A u] (g u] u 0 ) 2 + A a {g a a 0 ) 2 } 1 

k,F is the symmetric nuclear matter Fermi momentum 
and p = 2kp/3n 2 is the baryon number density. We have 
adjusted the symmetry energy to be e sym = 25.52 MeV 
at p = 0.1 fm - , and taking care that for the minimum 
possible values of L, the EOS of neutron matter has no 
binding. We point out that this value of the symmetry 
energy at 0.1 fm -3 is inside the range 25.5(10) MeV ob¬ 
tained in [38| from the properties of doubly magic nuclei. 


TABLE III: Values of some of the A i,g Pi pairs and their 
correspondent symmetry energy slopes. 


L (MeV) 



9p 

111 

A; 

0.00 

9.26 

75 

a ct 

0.01 

10.42 


Atj 

0.02 

10.18 

50 

a ct 

0.02 

12.76 



0.03 

11.78 


In Fig. [l] see also Table eh for the parameters A^, g Pi 
used, we plot the symmetry energy as a function of the 
density for three values of L, 111 MeV, corresponding to 
the model with no uip nor toa terms, 75 and 50 MeV. 
Below p = 0.1 fm , where all curves cross, the pcr term 
has a stronger effect, while above due to the saturation of 
the cr-field the pui term gives rise to a stronger softening. 
We expect, therefore, a stronger effect of the non-linear 
term pa on the properties of the crust while the pui term 
will have a larger influence on properties determined by 
the suprasaturation EOS. 

As reference, we plot in Fig. [2] the symmetry energy, 
the incompressibility K sym , the incompressibility coeffi¬ 
cient K t and the third derivative of the symmetry energy 
Q S ym , all defined at saturation, as a function of L. It is 
clearly seen the existence of correlations between the pa¬ 
rameters that define the isovector channel, and in general 
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these correlations agree with the correlations discussed 
in j27l. [28j. The non linear terms up and ap give similar 
results, although with the ap non-linear term a steeper 
behavior with L is obtained for both J and K T . 



FIG. 2: (Colour online) Symmetry energy and derivatives at 
saturation density as a function of the slope of the symmetry 
energy for the TM1-2 model. 


2. Hyperon couplings 

The NLWM extended lagrangian density includes also 
the mass of the eight baryons of the baryonic octet, the 
mass of a* and the (f> mesons, the electron and muon 
masses and the meson-hyperon couplings that we will 
discuss next. 

The u and <p isoscalar vector mesons coupling constants 
gu, B and are determined by SU(6) symmetry: 

111 . . 

t^9ojn 2 9^s 9^s (12a) 

2^/2 

‘2‘94'a ~ ^9<l>s ~ 9<t >s = ^ 9ojn- (l^b) 

For the p-meson we take 

9pn 9pt. 9pm 9 pa 0 ( 1 * 1 ) 

and the isospin operator in the baryon-meson coupling 
term in the Lagrangian takes into account the isospin 
symmetry. To define the hyperon-cr meson couplings g aB , 
we use the hyperonic potential in symmetric nuclear mat¬ 
ter at saturation defined by the relation 

U% = x u hVu - x aH V a (14) 

where V u = g u uo = 273.88 MeV, V a = g a ao = 
342.71 MeV are defined at saturation density, g aa = 
XaHg<j N and g UH = x^ng^N- While U f = -28 MeV 
is quite well constrained, Ui and U~ are not so well 
constrained: the experimental data [42j suggest that Ui 


is attractive, while Up is repulsive. In this work we es¬ 
tablished Ui = —28 MeV for all calculations and we 
use Ug = +30,0, -30 MeV, and Ug = +18,0, -18 MeV 
to take into account uncertainties. In Table m the op¬ 
tical potential sets (OPS) used in the present work are 
identified by numbers. 

TABLE IV: Reference table for the optical potentials sets used 
in the study. All sets satisfy Ui = —28 MeV. The optical 
potentials Up are defined in symmetric nuclear matter at 
saturation po- In the calculations which include the a * meson 
the values listed below for the strange a* meson couplings 
have been used. I 11 all cases we have = 8.524. In the last 
two lines the values of YY optical potentials Up at po /5 are 
given, where in all cases we have Fa(po/5) = —9.83 MeV. 
The potentials units are MeV. 


Set 

1 

2 

3 

4 

5 

6 

7 

8 

9 

up 

+30 

+30 

+30 

0 

0 

0 

-30 

-30 

-30 

Ui 

+18 

0 

-18 

+18 

0 

-18 

+18 

0 

-18 

Serf, 

9.87 

9.87 

9.87 

8.38 

8.38 

8.38 

6.10 

6.10 

6.10 

9crit 

13.01 

12.68 

12.27 

13.01 

12.68 

12.27 

13.01 

12.68 

12.27 

uE 

-15.75 

-15.74 

-15.74 

-8.80 

-8.80 

-8.80 

-1.81 

-1.81 

-1.81 

ui 

-9.57 

-7.37 

-5.24 

-9.57 

-7.37 

-5.24 

-9.57 

-7.37 

-5.24 


III. RESULTS AND DISCUSSION 

In the following we discuss the properties of the equa¬ 
tions of state and neutron stars obtained spanning the 
parameters presented in the last section. 

A. Equation of state of hyperonic stars 

In the present section we generalize the studies de¬ 
veloped in [20|, |37j, E3 and discuss the influence of the 
density dependence of the symmetry energy on the hy¬ 
peronic content of a neutron star. It has been shown 
in [2(1 [37] that star properties such as the radius, mass, 
strangeness content or the central baryonic density de¬ 
pend non-linearly on the slope L. 

To make a systematic study we cover the range 50 < 
L < 111 MeV, for both non-linear terms pu and pa, gen¬ 
erating an EOS and the corresponding family of neutron 
stars for each value of L. In Fig. [3] we show the stellar 
matter pressure at the saturation nuclear density, P(po), 
as a function of the slope L. The band, that restricts L 
to L < 88 MeV, identifies the allowed values of P(po) 
obtained in [44| from a microscopic calculation. In the 
following we will discuss the strangeness content in the 
star covering the whole range 50 < L < 111 MeV. 


We consider hyperon-hyperon interactions described 
by the exchange of the vector meson <f> and the scalar 
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FIG. 3: (Colour online) Stellar matter pressure at the satura¬ 
tion nuclear density po- The band defines the range of allowed 
values calculated in {44]] from a microscopic calculation. 

meson cr*, both with hidden strangeness. The A — A in¬ 
teraction was shown in [45| to be weakly attractive. This 
is implemented in RMF models by considering a small 
g a * B or even g a * B = 0 as described in Section IIIBI 

We have considered five different sets of coupling con¬ 
stants for the baryonic degrees of freedom in our nuclear 
matter EOS: a) one containing only neutrons, protons, 
electrons and muons in chemical equilibrium which is 
designated by np; b) in a second scenario, designated by 
no-(j>a* , hyperons are also included. The hyperon cou¬ 
pling constants to the vector mesons w and p are defined 
by Eqs. (112al) and (1131) . and the scalar meson a couplings 
are defined by Eq. (HI; c) a third scenario excludes the 
cr* but includes the exchange of the meson whose cou¬ 
pling constants are defined by the Eq. (112bl) with the 
SU(6) prescription. This case is designated by g<f, and 
takes into account that the recent results seem to indi¬ 
cate that the binding A — A is very weak; d) still keeping 
only the vector meson </> but lifting the SU(6) symmetry, 
the <^-hyperon couplings are defined by 

4^2 

= = <?<p= = ^ gun (15) 

This choice of couplings, designated by 2g c / > , brings ex¬ 
tra repulsion between hyperons, affecting mostly the 5- 
hyperon; e) a last possibility includes both the cr* and (f) 
mesons, keeping SU(6) symmetry to define the hyperon 
couplings to (j>. The g a * B couplings that have been used 
are listed in Table m This case is designated as cjxj*. 
The values chosen for g%, and listed in Table [TV] are such 
that taking for the $ the SU(6) values the respective sym¬ 
metric hyperon matter attains saturation for a binding 
energy of 1 MeV. The optical hyperon potentials Uy of 
the respective hyperon multiplet symmetric matter, 

u h = -gcrnVo - 9^0 + g UH u> 0 + (16) 

calculated at p = p 0 /5 are also fisted in IIVI 



0 1 2 3 4 5 

E (fm -4 ) 



1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

P/Po 



P/Po 

FIG. 4: (Colour online) Equation of state (top), strangeness 
fraction (middle) and electron fraction (bottom) for pn matter 
and different choices of the hyperon potentials for the TM1-2 
model. All figures have been done for U™ = +30 MeV and 
U~ = —18 MeV. The thick black lines correspond to L = 111 
MeV and the thin colored lines correspond to L = 50 MeV. 
Dots indicate the central energy density or baryonic density 
of the maximum mass star. 


The stellar matter EOS is obtained taking at low den¬ 
sities, below neutron drip, the Baym-Pethick-Sutherland 
[4(| EOS, which is adequate to describe the neutron star 
outer crust, between neutron drip and 0.01 fm -3 Bethe- 
Baym-Pethick EOS [4(| is considered, and above this 
density the EOS is smoothly interpolate to the homo¬ 
geneous equation of state which is taken above 0.5po- 
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TABLE V: Nucleonic maximum mass star properties. 



L 

M g 

M b 

R 

£0 

Pc 


111 

2.27 

2.67 

12.50 

5.216 

0.82 

A.oj 

50 

2.22 

2.63 

12.02 

5.474 

0.87 

a ct 

50 

2.24 

2.65 

12.09 

5.399 

0.86 


In Fig. []]we show the EOS, the strangeness fraction 
and the electron fraction for L = 50 and 111 MeV and the 
different choices of the couplings for the baryon degrees 
of freedom described above np, no-<f>cr*, g^, 2and 0cr*. 
The dots indicate the central energy or baryonic densities 
corresponding to the maximum mass neutron star. 

Some general conclusions, see also EMUS , can be 
drawn: a) both the strangeness and a small symmetry 
energy slope L give rise to a softer EOS. The inclusion 
of 0 excluding a* will produce a softer EOS than the 
pn-EOS but may still be quite hard if we do not restrict 
the choice of the couplings to the SU(6) symmetry; b) 
the inclusion of strangeness softens more an EOS with 
a large than with a small L because the fraction of hy- 
perons grows faster with density. For the same value of 
L, the up non-linear term favors more the appearance of 
hyperons than the top non-linear term; c) including only 
the 0-meson and excluding the a* reduces the fraction 
of hyperons due to the extra repulsion between hyperons 
introduced. If the a* is also included, the choice of the 
couplings will determine the behavior of the EOS and 
may soften the EOS more than when no strange meson 
is included; d) a larger amount of hyperons and a smaller 
L induce a smaller amount of electrons, and this effect 
is stronger for the EOS with the uip term. These main 
features will justify the results we discuss next for the 
dependence of the neutron star properties on the slope L 
and strangeness content. 


B. Dependence of hyperonic star properties on L 

The results obtained after integrating the TOV equa¬ 
tions, Eqs. (0), for all the EOS defined above are pre¬ 
sented and discussed in the present subsection. In Table 
m and Fig. 0 the maximum mass star properties ob¬ 
tained with an EOS excluding strangeness are given. The 
RMF model used allows maximum masses well above the 
limit imposed by the pulsar PSR J0348+0432 mass with 
M = 2.01 ± 0.04 Mg. Changing the symmetry energy 
slope, softens the EOS and as a result the central density 
is larger and the radius smaller. For the same value of L 
this effect is larger when the non-linear term up is used. 
Similar effect was discussed in filiH ]. The density depen¬ 
dence of the symmetry energy affects the maximum star 
mass only slightly: not more than 1% with the ap term, 
and < 3% with the uip term. In the last case the mini¬ 
mum mass is obtained for L ~ 80 MeV, while in the first 
the lowest value of L gives the smallest mass. The be- 



56 64 72 80 88 96 104 56 64 72 80 88 96 104 
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FIG. 5: (Colour online) Maximum gravitational mass (a), 
maximum mass star radius ( b ) versus L for np matter in j3 
equilibrium. 


havior obtained with the cop term is due to the interplay 
between a softening created for large L due to the larger 
amount of hyperons, and for small values of L due the 
smaller symmetry energy at large densities. 

In Fig. [6]we plot the families of stars obtained with the 
parametrizations 3 and 7 of Table IIVI and consider the 
different possibilities of including or not including 0 and 
a* mesons. For reference we also include the nucleonic 
stars. The dots mark the maximum mass configurations. 
The main conclusions drawn from the figure are effects 
already known: the presence of strangeness may have a 
strong effect on the maximum mass configuration, in par¬ 
ticular, a strong reduction occurs if the 0 meson is not 
included, or if it is included together with a* correspond¬ 
ing to a quite attractive YY potential, see, however, j47j 
where it was shown that if a strong enough g$ coupling 
is chosen, it was possible to describe a 2M© star includ¬ 
ing a* and 0 mesons; a small value of L will give rise to 
stars with a smaller radius, the effect being larger if the 
uip non-linear term is used. These general features are, 
however, sensitive to the choice of the hyperon optical 
potentials in symmetric nuclear matter as we will discuss 
next. Finally, strangeness is present in less massive stars 
if we take an attractive potential. 

Properties of maximum mass hyperonic stars are given, 
in appendix, in Tables IVIIIIVIIII and IIXI corresponding 
respectively to stars obtained with Uy, = +30, 0, —30 
MeV. Since there is a non-monotonic behavior of the 
maximum mass with L, in all cases we give information 
with respect to stars obtained with an EOS with L=110 
MeV, 50 MeV and the value of L that gives the star with 
the minimum maximum mass when this does not occur 
for L = 50 MeV. The mass Ms is the mass of a star 
with a droplet of strangeness in its center, i.e. it corre¬ 
sponds to the mass of the star that defines the onset of 
strangeness in the star core. 

In the two panels of Fig. [7] we plot, respectively, 













2.4 


2.12 



R (km) 

FIG. 6: (Colour online) Gravitational mass vs the radius for 
two different values of L with the parameters corresponding 
to set 3 (U£? = +30 MeV, top) and set 7 (U£ = —30 MeV, 
bottom). The dots indicate the maximum gravitational mass 
stars. The faded lines indicate the lower and upper limits of 
the pulsar PSR J0348+0432 mass. 

the gravitational mass (top) and the strangeness frac¬ 
tion (bottom) of the maximum mass star as a function 
of the slope L, taking U£ = +30 MeV. Results for hy- 
peronic stars are shown for the Ui? = +18,0,-18 MeV, 
the nonlinear terms up and uip, and different choices of 
the couplings to the <fi and a* mesons. 

All curves of plots (b) and (c) of the top panel fall 
within the mass limits (Mg = 1.97 ± 0.04 M 0 ) of the 
observed pulsar PSR J1614-2230 Q for any U~ with 
the nonlinear up term and almost all for the nonlin¬ 
ear ujp term, the exception being Uff = —18 MeV in 
plot c) in a small interval near the mass minimum. 
However, the mass of the pulsar PSR ,10348+0432 
Mg = 2.01 ± 0.04 Mq , is only attained with a stronger 
gif, than the one obtained imposing SU(6) symmetry. Ex¬ 
cluding the <j> and a* mesons, or choosing a too attrac¬ 
tive YY potential will not allow hyperonic stars to have 
a mass of the order of 2 Mg. The maximum masses ob¬ 




FIG. 7: (Colour online) Maximum gravitational mass and 
maximum mass star strangeness fraction versus L for Uj? = 
+30 MeV and all the other optical potentials (+18,0,-18) for 
Ui in MeV. 


tained with the attractive YY we have chosen are even 
smaller than the ones obtained excluding (f> and a*. We 
should point out that maximum mass configurations do 
not depend on US when the g^ is very strong because 
the H hyperon is disfavored due to its large strangeness 
charge and the repulsive effect of the </>-meson. 

Another conclusion that is drawn from the top panel 
of Fig. [7]is that the maximum mass has a non-monotonic 
behavior with L and generally it has a minimum for an 
intermediate value of L, just as was already seen for nu¬ 
cleonic stars in Fig. [5] A small L may give a larger max¬ 
imum mass due to a smaller strangeness content. This 
is the case when a* is excluded and g^ is defined by the 
SU(6) symmetry, see panels (c) and (g) in Fig. [TJ 

The total strangeness fraction evolves nonlinearly with 
L, see Fig. Q bottom panel, and may even show a non¬ 
monotonic behavior with L if no YY or an attractive YY 
interaction is chosen, figures a) and d). The large in¬ 
crease of the strangeness fraction for small values of L in 
bottom panel d) is due to the onset of the E° hyperon as 
is clearly seen in panels (d) and (h) of Fig. [8] obtained 
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FIG. 8: (Colour online) Onset particle density versus L with 
U£ = +30 MeV (both panels) and Ug = —18 MeV in the 
top panel and Ug = +18 MeV in the bottom panel. 


respectively for U~ = —18 and +18 MeV. An attractive 
U~ favors the onset of the 5“ hyperon which is the sec¬ 
ond hyperon to set in, see top panel. However, if g$ is 
too repulsive comes close with the E _ , see (b) top 
panel. A repulsive U£? together with a repulsive YY in¬ 
teraction will hinder completely the appearance of the 
S _ , see (f) and (g) of bottom panel. Fig. [8] also shows 
how the slope L influences the onset of the neutral versus 
negatively charged hyperons: the onset density of neutral 
hyperons decreases with increasing L , while an opposite 
behavior occurs for the negatively charges hyperons. A 
similar conclusion was drawn in [201]. 

As a general trend the larger the strangeness fraction 
the smaller the mass. This trend is broken if the g^ cou¬ 
pling is very strong and is also influenced by the behav¬ 
ior of the symmetry energy with the density. An asy-soft 
EOS will disfavor the strangeness onset if Uy is repulsive: 
in all plots of both panels of Fig. [51 the density of onset 
of the A-hyperon decreases when L increases. However, 
if the YY interaction is not too repulsive, intermediate 
values of L favor the appearance of negatively charged 
hyperons that decrease the electron fraction and soften 


the EOS, see Fig [4j Further decreasing L shifts the hy¬ 
peron onset to larger densities, which occurs with the A 
onset, and these effect is not compensated by the earlier 
onset of the negatively charged hyperons. 



L (MeVj L (MeV) 

FIG. 9: (Colour online) Maximum mass star radius versus 
L for TM1-2 with U™ = +30 MeV and all the other optical 
potentials (+18,0,-18) for Ug in MeV. 

In Fig. [9] the behavior of the maximum mass star 
radius with L is compared for the different hyperon 
parametrizations considered taking = +30 MeV. The 
general trend is a decrease of the radius if L decreases. 
This radius reduction may be as large as 1 to 1.5 km if 
no YY or a quite attractive YY interaction is considered. 

For the parametrization labelled cjxr* the hyperons A, 
E _ and are present in the core independently of L. 
Let us now consider U~ = +18 MeV, which gives the 
smallest radii and has the largest fractions of strangeness 
with the nonlinear term up. For L ~ 90 — 95 MeV the 
curves for the radius and strangeness fraction suffer a 
faster change, respectively increase and decrease, when 
L increases. This is due to the onset of the E° hyperon 
and occurs for both up and crp nonlinear interactions, 
see Fig. [5] A similar situation occurs for Ui? = 0 and 
— 18 MeV, however in this case E° sets in only for the up 
mixture and for L < 80 MeV. The radii and strangeness 
fraction in this case come close to the ones obtained with 
US = +18 MeV. It is interesting to notice that this is 
a quite different behavior from the trend obtained with 
the up nonlinear term where the radius difference be¬ 
tween L = 50 and 110 MeV is not larger than 1 km and 
the strangeness fraction for L = 50 MeV is half the one 
obtained with the up term, see plot (h) of the bottom 
panel of Fig. [7] and plot (d) of Fig. [9] 

In Fig. [10] the strangeness fraction is plotted for Uj? = 
—30 MeV. In this case the E _ hyperon is the first to set 
in. Since a smaller L favors the formation of negatively 
charged hyperons, the general trend is a decrease of the 
strangeness when L increases, even for a repulsive YY 
interaction. For the attractive YY parametrization the 
fast decrease of the strangeness for L > 90 MeV is due 
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FIG. 10: (Colour online) Maximum mass star strangeness 
fraction versus L for TM1-2 with = —30 MeV and all the 
other optical potentials (+18,0,-18) for US in MeV. 


to the disappearance of the £° hyperon. 


FIG. 12: (Colour online) Star mass of onset of strangeness 
versus L for TM1-2 for U~ = —18 MeV and several values 
for U™. For U£ = 0, the full dots indicate the critical slope, 
L = 72 MeV for A w and L = 67.2 MeV for A CT , above which 
Pe onset > PAonset • For Ug = -30 (+30) the first hyperon to 
appear is the E“. 



L (MeV) 


FIG. 11: (Colour online) Competition to first hyperon onset 
between the E _ and the A elements of the lightest baryon 
octet when we change L+ . 

Fig. [TT1 clarifies the existing competition between the 
A and the £“ onset. An attractive U£ potential favors 
the £“ onset. On the contrary a repulsive favors 
the A, but for Uj? = 0 or close the L defines the hy¬ 
peron onset: small L favors the £“. This behavior was 
already discussed in [20]. The other parametrizations of 
the hyperon-meson couplings give rise to a similar behav¬ 
ior and, therefore, are not represented. 

The star mass of strangeness onset, i.e., the mass of 
the star with a droplet of strangeness at the center and 
an infinitesimal Ns/Nb, is plotted in Fig. [T3] 

These curves do not depend on the inclusion of strange 
mesons, which only act in matter with a finite amount 
of hyperons. The A and £ _ hyperons compete to ap¬ 
pear first as and L change. The repulsive = +30 
MeV or even the U£ = 0 MeV potential favor the on¬ 


set of the A hyperon, as already referred above. How¬ 
ever, in the last case, for a low enough L, the £~ sets 
in first if L < 67.2 MeV for ap or L < 72 MeV for the 
cop non-linear term. Below those critical values of L the 
strangeness onset mass decreases with the decrease of L. 
If = —30 MeV the onset of the £~ occurs always at 
a smaller density than A. If the first hyperon to set in 
is the A the star mass of strangeness onset is ~ 1.4 Mq 
and the dependence on the slope L is small. However, if 
is attractive, £ - may be the first hyperon to set in 
first, and, in this case, the strangeness onset star mass is 
sensitive to L. For = —30 MeV this mass decreases 
from ~ 1.3 to 1.05 Mq when L decreases from 111 to 
50 MeV. We conclude that if the potential is repul¬ 
sive as experiments seem to indicate [42| we may expect 
strangeness in stars that have a mass at least as large as 
~ 1.4Mq. Less massive stars are totally determined by 
the nucleonic properties of the EOS. 


We have seen that strangeness is present in stars with 
masses (above or equal) to ~ 1M 0 — 1.42M 0 depend¬ 
ing on . L strongly influences neutron star radius 
and the strangeness content. To understand the effect 
of strangeness and L on the star radius, we have fixed 
the star mass and calculated the radius changing both L 
and strangeness. In Fig. [13] we plot the results obtained 
for a 1.67 M 0 and a 1.75 M 0 star. It is clearly seen 
that the radius decreases linearly with the increase of 
the strangeness fraction for all the nucleonic EOS shown, 
which differ in the value of L and non-linear p term in the 
Lagrangian. We have fitted the data for the star radius 
as a function of the strangeness fraction through 

R = rjsY s + Ro (km), 


(17) 
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FIG. 13: (Colour online) Radii versus total strangeness con¬ 
tent of the stars with masses 1.67Mq and 1.75Mq for L = 111 
MeV (circles, full line), L = 75 MeV (triangles), L = 50 MeV 
(squares). The straight lines are fitted to the data obtained. 
The red symbols and red dotted line refers to the EOS with 
the ap term and the blue symbols and dashed-dotted to the 
EOS with the cup term. The points represented include infor¬ 
mation from all the hyperon-meson couplings considered. 


TABLE VI: Parameters (slope r/s = d ^ N ^, I J' NB ' ) an d zero 
strangeness radius Ro = R(Ys = 0)) for the straight lines that 
describe the star radius versus the total strangeness fraction 
for four values of the star mass: 1.44, 1.60, 1.65 and 1.75 Mq, 
see Fig. 1131 The uncertainties in the fits are below 6% in the 


slope and 

below 

0.02% 

in R 0 







1.44 

Mq 

1.60 

Mq 

1.67 

Mq 

1.75 

Mq 

L 

Vs 

Ro 

Vs 

Ro 

Vs 

Ro 

Vs 

Ro 

(MeV) 

(km) 

(km) 

(km) 

(km) 

(km) 

(km) 

(km) 

(km) 

HI 

10.30 

14.32 

10.69 

14.23 

10.72 

14.19 

10.22 

14.11 

75 (A ct ) 

10.82 

13.66 

10.09 

13.63 

10.46 

13.61 

10.47 

13.58 

75 (A„) 

10.84 

13.59 

10.20 

13.55 

10.21 

13.52 

10.87 

13.50 

50 (A ct ) 

10.22 

13.32 

10.16 

13.34 

10.12 

13.34 

10.26 

13.32 

50 (A„) 

10.33 

13.27 

10.31 

13.28 

10.13 

13.26 

10.17 

13.24 


where R 0 = R(Y S = 0) and rjs = Y s = N S /N B - 
The parameters of the fitted lines for four values of the 
star mass are listed in Table eh The values of the slope 
are approximately equal within the uncertainty of the 
slope and agree with the values obtained previously in 

m 

Radius estimates inferred from photospheric radius ex¬ 
pansion bursts and thermal emissions from quiescent low- 
mass X-ray binaries and isolated neutron stars using 
Bayesian techniques [H] indicate that the radius of a 
1.4Mq star is 12.1±1.1 km. If 1.4Mq stars are purely 
nucleonic the radius gives information on the density de¬ 
pendence of the symmetry energy, and small values of L 
are favored. According to our study if the U£ is repul¬ 
sive, it is very likely that 1.4M 0 stars are nucleonic stars. 
However, if is attractive, small values of L favor an 


early appearance of strangeness and the radius is deter¬ 
mined by both the density dependence of the symmetry 
energy and the strangeness content. 


IV. CONCLUSIONS 

In the present study we discuss the joint effect of the 
density dependence of the symmetry energy and hyperon 
interaction on several properties of neutron stars, includ¬ 
ing radius, gravitational mass, strangeness content and 
star mass at strangeness onset. This was carried out 
within the relativistic mean field framework. The den¬ 
sity dependence of the symmetry energy was modelized 
including a non-linear up or ap term in the Lagrangian 
density. We have changed the symmetry energy slope 
at saturation between 50 and 110 MeV, taking care that 
the parametrizations obtained do no predict neutron star 
matter that saturates at a finite density. 

For the hyperon meson interaction we have considered 
the isoscalar vector mesons w,^-coupling constants de¬ 
fined by the SU(6) symmetry, while the p-meson cou¬ 
pling is defined by the isospin symmetry of the hyperon, 
taking the nucleon coupling as reference. In order to 
get a quite repulsive YY interaction, we have also con¬ 
sidered a parametrization with broken SU(6) symmetry, 
and took a parametrization with g$ > g l p(SU(6)). To 
define the hyperon-er meson couplings we use the hyper- 
onic potential in symmetric nuclear matter at saturation 
and take Uj^ = —28 MeV. Since and U§ are not 
well constrained we have considered several possibilities 
from attractive to repulsive potentials for these two hy- 
peronic potentials. The hyperon-cr* coupling was chosen 
non-zero in a single parametrization which gives rise to 
a quite attractive YY interaction. 

Next we summarize some of the main conclusions. 

Nucleonic stars have smaller radii for EOS with smaller 
values of L as shown in [32|, and this effect is larger if 
the non-linear up instead of the ap term is used. These 
two terms implement different density dependencies on 
the symmetry energy due to the different behavior of the 
a and u mesons, the first one saturates at high densities 
while the second one increases almost linearly with den¬ 
sity. Consequently the ap (up) term has a stronger effect 
below (above) saturation density. The gravitational mass 
is not affected by more than ~ 3% when L changes from 
50 to 110 MeV. 

Including hyperons will soften the EOS, however 
quantitatively the effect depends a lot on the hyperon 
parametrization chosen. We have shown that 2 Mq are 
obtained if the YY interaction is repulsive enough, as 
shown in previous works p, D3. GH 0, S3 ■ Within the 
TM1 parametrization chosen, it could be shown that if 
the hyperonic potential is repulsive there will only 
exist strangeness in the core for stars with M > 1.4Mq, 
and that the first hyperon to set in is the A. On the other 
hand an attractive U£ favors the onset of the E _ before 
the A and a 1M 0 star may already have strangeness in 











12 


its core. If the potential U^f is close to zero the value of 
L will define whether the A or the £“ will first set in, 
with a small L favoring the appearance of the £ _ . 

We could confirm the results of [2(| that star radius 
depends linearly on the total strangeness fraction, and 
we have shown that the slope of this linear correlation 
does not depend neither in L nor in the star mass and is 
~ 10.5 knr/ 775 . 

The existence of strangeness in the core of a 1.4 Mq 
star may have to be considered when the £ potential 
in nuclear matter is not repulsive. In this case, the as- 
trophysical observations will give information not only 
on the slope L of the symmetry energy but also on the 
strangeness content that may be estimated using expres¬ 
sion (fT71) . 

We have also shown that while in nucleonic stars re¬ 
ducing the slope L from 110 to 50 MeV may give rise to 
a decrease of the maximum mass star radius of ~ 0.5 km, 
for hyperonic stars this reduction can be of the order of 
1 km or even larger if an attractive YY is considered. 

Finally it should also be referred that to understand 
the structure of a neutron star it is not enough to con¬ 
strain properties at saturation density, but information 
on the density dependence of the equation of state, in par¬ 
ticular, the symmetry energy is also needed. Taking two 
different parametrizations for the isospin channel of the 


equation of state it was shown the density dependence of 
the symmetry at high density influences properties such 
as the mass, radius or strangeness content of the star. 

In this study we were very conservative and we have 
only considered for the hyperon couplings one possibility 
that did not satisfy the SU( 6 ) symmetry. In f47j several 
parametrizations have been proposed that do not satisfy 
either the SU( 6 ) symmetry and even with a YY attractive 
interaction in some of the channels it was possible to 
describe a 2 Mq star. With the present work we just want 
to show in a more systematic way that it is important to 
put constraints on the hyperon interactions before the 
interpretation of observational data may be used to set 
constraints on the equation of state. 
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Appendix 

In this section we display the properties of maximum 
mass stars calculatd with the EOS presented. 
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TABLE VII: Maximum star properties for the TM1-2 model. The hyperon potentials are U£ = +30 MeV and U£ = 
+ 18, 0, —18 MeV. The minimum gravitational mass star is also listed in the intermediate critical L c . Ms is the star mass of 
onset of strangeness for a certain L. 



L 

M g 

U£ = +18 
M b R £ 0 

Vs 

Pc 

M s 
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M g 

M b 

II 
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Vs 
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L 

M g 

M b 

II 

-18 

£o 

Vs 

Pc 

Ms 












no-(j>a* 














111 

1.80 

2.02 

13.19 

4.30 

0.087 

0.75 

1.41 


1.77 

1.99 

13.34 

4.07 

0.077 

0.72 

1.41 


1.72 

1.92 

13.43 

3.92 

0.076 

0.69 

1.41 

A co 

96 

1.79 

2.01 

12.88 

4.51 

0.088 

0.78 

1.38 


1.75 

1.97 

12.99 

4.34 

0.083 

0.76 

1.37 


1.70 

1.91 

13.00 

4.28 

0.091 

0.75 

1.38 


50 

1.81 

2.07 

12.47 

4.67 

0.079 

0.81 

1.41 


1.79 

2.03 

12.56 

4.49 

0.077 

0.78 

1.41 


1.74 

1.96 

12.50 

4.535 

0.093 

0.79 

1.41 

K 

96 

1.79 

2.02 

12.95 

4.45 

0.088 

0.77 

1.38 


1.76 

1.98 

13.06 

4.26 

0.082 

0.75 

1.38 


1.71 

1.92 

13.09 

4.18 

0.088 

0.74 

1.38 


50 

1.82 

2.07 

12.55 

4.57 

0.079 

0.79 

1.41 


1.79 

2.03 

12.64 

4.40 

0.076 

0.77 

1.41 


1.74 

1.97 

12.60 

4.41 

0.090 

0.77 

1.41 


111 

2.12 

2.45 

12.31 

5.48 

0.078 

0.87 

1.39 


2.12 

2.45 

^9<t> 

12.30 

5.48 

0.078 

0.87 

1.42 


2.12 

2.45 

12.30 

5.48 

0.078 

0.87 

1.39 

A (jj 

79 

2.08 

2.42 

12.01 

5.63 

0.067 

0.90 

1.38 


2.08 

2.42 

12.01 

5.63 

0.067 

0.90 

1.38 


2.08 

2.42 

12.01 

5.63 

0.067 

0.90 

1.38 


50 

2.09 

2.44 

11.91 

5.60 

0.062 

0.90 

1.41 


2.09 

2.44 

11.91 

5.60 

0.062 

0.90 

1.41 


2.09 

2.44 

11.92 

5.60 

0.062 

0.90 

1.41 

K 

68 

2.10 

2.45 

12.03 

5.57 

0.068 

0.89 

1.38 


2.10 

2.45 

12.03 

5.57 

0.068 

0.89 

1.38 


2.10 

2.45 

12.03 

5.57 

0.068 

0.89 

1.38 


50 

2.11 

2.46 

11.97 

5.54 

0.064 

0.89 

1.41 


2.11 

2.46 

11.97 

5.54 

0.064 

0.89 

1.41 


2.11 

2.49 

11.97 

5.54 

0.065 

0.89 

1.41 


111 

1.80 

2.02 

13.19 

4.30 

0.087 

0.75 

1.41 


1.96 

2.23 

94> 

12.44 

5.31 

0.125 

0.87 

1.41 


1.94 

2.21 

12.41 

5.36 

0.138 

0.88 

1.41 

A ijj 

86 

1.94 

2.22 

12.14 

5.50 

0.111 

0.90 

1.37 


1.94 

2.22 

12.14 

5.48 

0.113 

0.90 

1.37 

90 

1.92 

2.19 

12.10 

5.58 

0.133 

0.92 

1.36 


50 

1.81 

2.07 

12.47 

4.67 

0.079 

0.81 

1.41 


1.95 

2.25 

11.99 

5.44 

0.100 

0.90 

1.41 


1.94 

2.23 

11.94 

5.53 

0.118 

0.91 

1.40 

K 

80 

1.96 

2.24 

12.18 

5.44 

0.112 

0.90 

1.38 

89 

1.95 

2.23 

12.23 

5.41 

0.119 

0.89 

1.38 

92 

1.93 

2.21 

12.18 

5.51 

0.137 

0.91 

1.38 


50 

1.82 

2.07 

12.55 

4.57 

0.079 

0.79 

1.42 


1.96 

2.26 

12.06 

5.38 

0.103 

0.89 

1.41 


1.95 

2.24 

11.99 

5.47 

0.121 

0.90 

1.41 


111 

1.80 

2.02 

13.19 

4.30 

0.087 

0.75 

1.41 


1.61 

1.78 

13.85 

3.30 

0.038 

0.60 

1.39 


1.58 

1.74 

13.94 

3.15 

0.032 

0.58 

1.39 

A U 

89 

1.60 

1.78 

12.78 

4.40 

0.108 

0.78 

1.38 


1.58 

1.76 

13.27 

3.66 

0.052 

0.66 

1.37 


1.55 

1.72 

13.32 

3.58 

0.053 

0.65 

1.35 


50 

1.81 

2.07 

12.47 

4.67 

0.079 

0.81 

1.41 


1.61 

1.80 

12.31 

4.68 

0.119 

0.83 

1.41 


1.59 

1.77 

12.28 

4.73 

0.130 

0.83 

1.40 

K 

89 

1.61 

1.79 

12.98 

4.10 

0.090 

0.73 

1.39 


1.59 

1.77 

13.39 

3.56 

0.048 

0.65 

1.38 


1.56 

1.73 

13.46 

3.45 

0.047 

0.63 

1.36 


50 

1.82 

2.07 

12.55 

4.57 

0.079 

0.79 

1.39 


1.62 

1.81 

12.86 

3.82 

0.058 

0.69 

1.39 


1.59 

1.78 

12.84 

3.81 

0.064 

0.69 

1.41 
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TABLE VIII: Maximum star properties for the TM1-2 model. The hyperon potentials are Ui = 0 MeV and US = +18, 0, —18 
MeV. The minimum gravitational mass star is also listed in the intermediate critical L c . Ms is the star mass of onset of 
strangeness for a certain L. 





Ui = +18 






ui = 

= 0 






ui = 

-18 





L 

M g 

M b 

R 

£o 

Vs 

Pc 

M s L 

M g 

M b 

R 

£o 

Vs 

Pc 

Ms 

L 

M g 

M b 

R 

£o 

Vs 

Pc 

M s 











no-(j>a* 














111 

1.78 

2.00 

13.06 

4.44 

0.098 

0.77 

1.41 

1.76 

1.98 

13.27 

4.15 

0.081 

0.73 

1.41 


1.72 

1.92 

13.42 

3.93 

0.076 

0.70 

1.41 


86 

1.74 

1.96 

12.43 

5.02 

0.118 

0.86 

1.37 

1.73 

1.95 

12.60 

4.75 

0.103 

0.82 

1.37 

90 

1.70 

1.90 

12.82 

4.48 

0.099 

0.78 

1.37 


50 

1.76 

1.99 

12.09 

5.19 

0.119 

0.89 

1.35 

1.75 

1.99 

12.22 

4.96 

0.107 

0.86 

1.35 


1.73 

1.95 

12.33 

4.78 

0.103 

0.83 

1.35 

a ct 

83 

1.75 

1.98 

12.51 

4.89 

0.115 

0.84 

1.37 86 

1.74 

1.96 

12.74 

4.59 

0.099 

0.80 

1.38 

92 

1.71 

1.91 

12.99 

4.28 

0.092 

0.75 

1.38 


50 

1.77 

2.00 

12.19 

5.05 

0.116 

0.87 

1.38 

1.76 

1.99 

12.33 

2 9<t> 

12.31 

4.841 

0.104 

0.84 

1.38 


1.73 

1.96 

12.46 

4.62 

0.099 

0.81 

1.38 


111 

2.11 

2.44 

12.31 

5.47 

0.080 

0.87 

1.41 

2.11 

2.44 

5.47 

0.080 

0.87 

1.41 


2.11 

2.44 

12.31 

5.47 

0.080 

0.87 

1.41 

A oj 

73 

2.05 

2.38 

11.88 

5.78 

0.076 

0.93 

1.37 

2.05 

2.38 

11.88 

5.78 

0.076 

0.93 

1.38 


2.05 

2.38 

11.89 

5.78 

0.076 

0.93 

1.38 


50 

2.06 

2.40 

11.79 

5.77 

0.074 

0.93 

1.35 

2.06 

2.40 

11.79 

5.77 

0.074 

0.93 

1.35 


2.06 

2.40 

11.79 

5.77 

0.074 

0.93 

1.35 

a ct 

50 

2.08 

2.42 

11.86 

5.69 

0.076 

0.92 

1.38 

2.08 

2.42 

11.86 

5.69 

0.076 

0.92 

1.38 


2.08 

2.42 

11.86 

5.69 

0.076 

0.92 

1.38 


111 

1.94 

2.21 

12.40 

5.37 

0.129 

0.88 

1.41 

1.94 

2.21 

9<t> 

12.40 

5.37 

0.128 

0.88 

1.41 


1.93 

2.20 

12.42 

5.35 

0.134 

0.88 

1.40 

A„ 

79 

1.90 

2.17 

11.92 

5.76 

0.127 

0.95 

1.37 

1.90 

2.17 

11.92 

5.76 

0.127 

0.95 

1.37 

84 

1.90 

2.16 

11.96 

5.73 

0.129 

0.94 

1.37 


50 

1.91 

2.19 

11.78 

5.77 

0.121 

0.95 

1.33 

1.91 

2.19 

11.77 

5.77 

0.121 

0.95 

1.33 


1.91 

2.19 

11.78 

5.76 

0.121 

0.95 

1.33 

a ct 

71 

1.91 

2.19 

11.97 

5.66 

0.127 

0.93 

1.39 

1.91 

2.19 

11.96 

5.66 

0.127 

0.93 

1.39 

75 

1.91 

2.19 

12.01 

5.64 

0.130 

0.93 

1.38 


50 

1.92 

2.20 

11.85 

5.67 

0.123 

0.93 

1.36 

1.92 

2.20 

11.85 

<j)a * 

13.79 

5.67 

0.123 

0.93 

1.3 


1.92 

2.20 

11.85 

5.66 

0.123 

0.93 

1.36 


111 

1.59 

1.76 

13.44 

3.78 

0.079 

0.68 

1.39 

1.59 

1.76 

3.33 

0.044 

0.61 

1.40 


1.57 

1.73 

13.98 

3.09 

0.027 

0.57 

1.39 

A„, 

91 

1.56 

1.73 

12.39 

5.03 

0.154 

0.88 

1.35 

1.55 

1.72 

12.61 

4.68 

0.132 

0.83 

1.37 


1.53 

1.69 

12.81 

4.34 

0.115 

0.77 

1.35 


50 

1.59 

1.78 

11.83 

5.54 

0.171 

0.96 

1.34 

1.58 

1.77 

12.00 

5.24 

0.152 

0.91 

1.33 


1.56 

1.74 

12.11 

5.02 

0.144 

0.88 

1.33 

a ct 

89 

1.57 

1.74 

12.55 

4.79 

0.141 

0.84 

1.38 

1.56 

1.73 

12.78 

4.42 

0.118 

0.79 

1.38 


1.54 

1.70 

13.03 

4.02 

0.094 

0.72 

1.38 


50 

1.59 

1.78 

11.95 

5.32 

0.163 

0.92 

1.36 

1.59 

1.78 

11.95 

5.32 

0.163 

0.92 

1.36 


1.57 

1.75 

12.24 

4.81 

0.135 

0.85 

1.35 
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TABLE IX: Maximum star properties for the TM1-2 model. The hyperon potentials are U£ = —30 MeV and US = +18, 0, —18 
MeV. The minimum gravitational mass star is also listed in the intermediate critical L c . Ms is the star mass of onset of 
strangeness for a certain L. 



L 

M g 

u£ = - 

M b R 

-f 18 

£o 

Vs 

Pc 

M s L 

M g 

M b 

II 

I ? 05 

- 0 

£o 

Vs 

Pc 

Ms 

L 

M g 

U£ = 
M b R 

-18 

£o 

Vs 

Pc 

M s 











nc 

i ~(j>a* 














111 

1.69 

1.89 

13.11 

4.37 

0.101 

0.77 

1.29 

1.69 

1.89 

13.11 

4.37 

0.101 

0.77 

1.29 


1.68 

1.87 

13.21 

4.23 

0.093 

0.74 

1.30 


86 

1.65 

1.84 

12.28 

5.22 

0.140 

0.90 

1.13 

1.65 

1.84 

12.28 

5.22 

0.140 

0.90 

1.13 

89 

1.64 

1.83 

12.41 

5.07 

0.133 

0.88 

1.14 


50 

1.67 

1.88 

11.88 

5.49 

0.146 

0.94 

1.02 

1.67 

1.88 

11.87 

5.49 

0.146 

0.94 

1.02 


1.67 

1.88 

11.91 

5.47 

0.144 

0.94 

0.99 

a ct 

83 

1.66 

1.86 

12.38 

5.06 

0.136 

0.87 

1.13 

1.66 

1.86 

12.39 

5.06 

0.136 

0.87 

1.13 


1.65 

1.85 

12.49 

4.93 

0.129 

0.85 

1.14 


50 

1.67 

1.89 

11.97 

5.37 

0.145 

0.92 

1.04 

1.67 

1.89 11.97 

5.37 

0.145 

0.92 

1.0 


1.67 

1.88 

12.01 

5.30 

0.141 

0.91 

1.02 


111 

2.09 

2.41 

12.16 

5.65 

0.096 

0.90 

1.30 

2.09 

A 9<t> 

2.41 12.16 

5.65 

0.096 

0.90 

1.28 


2.09 

2.41 

12.16 

5.65 

0.096 

0.90 

1.30 

Ac j 

70 

2.02 

2.33 

11.61 

6.13 

0.102 

0.98 

1.05 

2.02 

2.34 

11.59 

6.14 

0.102 

0.98 

1.04 


2.02 

2.33 

11.60 

6.13 

0.102 

0.98 

1.05 


50 

2.02 

2.35 

11.51 

6.15 

0.102 

0.98 

0.99 

2.02 

2.35 

11.51 

6.15 

0.102 

0.98 

0.99 


2.02 

2.35 

11.51 

6.15 

0.102 

0.98 

0.99 

a ct 

50 

2.04 

2.37 

11.59 

6.05 

0.104 

0.97 

1.01 

2.04 

2.37 

11.59 

6.05 

0.104 

0.97 

1.01 


2.04 

2.37 

11.58 

6.05 

0.104 

0.97 

1.01 


111 

1.89 

2.15 

12.24 

5.584 

0.152 

0.92 

1.30 

1.89 

2.15 

9<t> 

12.24 

5.58 

0.152 

0.92 

1.30 


1.89 

2.15 

12.24 

5.58 

0.152 

0.92 

1.30 

A„ 

84 

1.84 

2.09 

11.60 

6.252 

0.166 

1.02 

1.11 

1.84 

2.09 

11.60 

6.25 

0.166 

1.02 

1.11 


1.84 

2.09 

11.60 

6.25 

0.166 

1.02 

1.11 


50 

1.85 

2.12 

11.38 

6.332 

0.163 

1.03 

1.01 

1.85 

2.12 

11.39 

6.33 

0.163 

1.03 

1.01 


1.85 

2.12 

11.38 

6.33 

0.163 

1.03 

1.01 

a ct 

73 

1.86 

2.12 

11.63 

6.150 

0.167 

1.00 

1.08 

1.86 

2.12 

11.62 

6.15 

0.167 

1.00 

1.08 


1.86 

2.12 

11.63 

6.15 

0.167 

1.00 

1.08 


50 

1.86 

2.13 

11.47 

6.207 

0.165 

1.01 

1.01 

1.86 

2.13 11.48 

6.21 

0.165 

1.01 

1.04 


1.86 

2.13 

11.47 

6.21 

0.165 

1.01 

1.01 


111 

1.57 

1.74 

12.53 

5.035 

0.180 

0.88 

1.27 

1.56 

1.73 13.06 

4.27 

0.122 

0.76 

1.27 


1.54 

1.70 

13.52 

3.65 

0.075 

0.66 

1.29 

A„ 

102 

1.56 

1.73 

11.65 

6.375 

0.261 

1.07 

1.21 

1.55 

1.71 

12.02 

5.70 

0.219 

0.98 

1.19 


1.52 

1.68 

12.40 

5.07 

0.181 

0.89 

1.19 


50 

1.60 

1.79 

10.95 

7.100 

0.272 

1.18 

0.99 

1.59 

1.78 

11.17 

6.64 

0.249 

1.11 

1.00 


1.57 

1.76 

11.30 

6.38 

0.241 

1.08 

1.00 

a ct 

102 

1.57 

1.73 

11.88 

5.982 

0.241 

1.02 

1.23 

1.55 

1.72 

12.32 

5.23 

0.189 

0.91 

1.23 

99 

1.53 

1.69 

12.60 

4.86 

0.166 

0.85 

1.21 


50 

1.60 

1.79 

11.04 

6.930 

0.270 

1.15 

1.03 

1.59 

1.78 

11.29 

6.42 

0.243 

1.08 

1.04 


1.57 

1.75 

11.45 

6.16 

0.231 

1.04 

1.01 





